Rectification in single molecular dimers with strong polaron effect 
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We study theoretically the transport properties of a molecular two level system with large electron- 
vibron coupling in the Coulomb blockade regime. We show that when the electron- vibron coupling 
induces polaron states, the current-voltage characteristic becomes strongly asymmetric because, in 
one current direction, one of the polaron state blocks the current through the other. This situation 
occurs when the coupling between the polaron states is smaller than the coupling to the leads. We 
discuss the relevance of our calculation for experiments on C140 molecules. 



I. INTRODUCTION 

The possibility of designing molecular junctions that 
behave as rectifiers was suggested by Aviram and Ratner 1 
based on an asymmetric two-level donor molecular sys- 
tem, with one level empty and one filled (an acceptor 
and a donor level). In one bias direction the levels 
are tuned closer to resonance, whereas for the reverse 
bias they are further detuned, giving rise to an asym- 
metric current-voltage (IV) characteristic. This scenario 
has been investigated in molecular systems for a num- 
ber of years^ The transport properties depend crucially 
on both the Coulomb interactions and the relaxation of 
the nuclei during the electron transfer processes. The 
Coulomb interactions limit the allowed charge configura- 
tions, while the electron- vibron coupling determines the 
overlap between the different configurations as well as 
provides channels for energy relaxation. 

Molecular transistors with strong electron- 
vibron coupling has recently become an ac- 
tive research area, both experimentall^ 3 ^^ and 
theoretically 7 i 8 i 9 i 10 i n i 12 i 13 i 14 i 15 i 16 i 17 Different types of 
vibrational modes have been observed. In the original 
work of Park et al^ clear signatures of phonon sidebands 
were seen in the tunnel spectrum. The size of the indi- 
vidual steps in the IV curves agreed well with a simple 
model based on Franck-Condon physics. 3 Furthermore, 
sidebands caused by internal vibration of the molecular 
systems have also been identified experimentally and, 
in particular, in a recent experiment by Pasupathy et 
alA where tunneling through dimerized C70 molecules 
was measured. These authors claimed that a vibrational 
mode at 11 meV due to the relative motion of the two 
C70 molecules was seen in the experiment. 

In this paper, we study theoretically transport through 
a dimer molecule with a single internal vibration. As a 
generic model, we include a coupling between the charge 
on the two parts of the dimer. The electron- vibron cou- 
pling thus introduces a possibility of forming polaronic 
states. Transport in systems with polaronsiS is a well- 
studied subject and the polaron formation is known to 
reduce the mobility. In single-electron transistor systems 
polaronic effects have also been studied theoretically, in 
particular for single- level S y S tems 7 *2i2ii2iAAii2iiiiiii£ , and 
only recently for a two level system^ Here we study a 
similar system but with the important modification that 
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FIG. 1: Dimer model studied in this paper. The two parts 
of the dimer are represented by single electronic levels, a and 
b, connected by a tunneling amplitude, t. The intramolecular 
vibration is illustrated as a spring connecting the two levels. 
A force acts on the part of the dimer that is occupied by an 
electron due to an electric field in the gap. 



the two levels in addition are coupled by tunneling, in 
which case an internal polaron can form. The main point 
of our study is that when this electron- vibron coupling is 
sufficiently strong and a polaron state is formed, it may 
lead to rectification effects and negative differential con- 
ductance, depending on how the molecule is situated in 
the constriction. 

The paper is organized as follows. In Section II we 
set up the model describing the dimer and the electron- 
vibron coupling, in Section III we discuss the physical 
concepts in a semi-classical picture, in Section IV we cal- 
culate the transport in the situation where the tunnel- 
ing rates is smaller than the polarons couplings, whereas 
Section V discusses the situation when the polarons are 
localized on the time scale set by tunneling to the leads. 
The cross-over between these two regime is studied in 
Section VI and finally Section VII discusses the experi- 
mental relevance and concludes. 



II. THE MOLECULAR DIMER MODEL 

The molecular model that we study (see Fig. is a 
dimer with 2 electron sites, or in a different terminology 
a "double quantum dot" . For simplicity, we neglect both 
spin and possible orbital degeneracies in the model. The 
two sites are coupled by tunneling. Furthermore, due 
to an electric field across the gap the internal molecular 
vibrational mode couples to the charges difference on the 
two sites. The electric field is produced either by the bias 
voltage, local charge traps, or image charges created in 
the leads by the charged molecule. In total our model 
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Hamiltonian for the molecule reads 
p 2 1 

H ° = 2^ + 2 mCj2 ^ 2 + A ^ a n ^ + + 

+ ^ ™a™6 - ^(4 C 6 + C l C a) + A Oa ™bK (!) 

Here ct (c^) (j=a,b) are the electron creation (anni- 
hilation) operators, whereas b^(b) is the vibron cre- 
ation(annihilation) operator and rij = CjCy Further- 
more, V g is the gate voltage which couples to the total 
charge on the molecule, A describes the energy difference 
between the two sites and, finally, U the Coulomb repul- 
sion between electrons occupying the two levels. The 
total Hamiltonian describing the molecule, the leads and 
the coupling between them is 

H = H D + (H n + H Tv ), (2) 

rj=LM 

where 

Hf] ~ ^ ^krj c krj C krj^ V ~ L,R, (3) 
k 

H Tr) = (^kr),j C kr] C j + T krj,j C j C fery) • ( 4 ) 

k,j=a,b 

Throughout, we will assume large U and the gate volt- 
age tuned such that the total occupancy is either zero or 
one. The subspace with zero electrons is, of course, easily 
diagonalized. With one electron on the dimer, the situa- 
tion is more complicated. Introducing Pauli matrices for 
the electronic degree of freedom and the dimensionless 
parameters (using h = 1) 

A2 A2 a A m 

<* = — 27, g = - — o, o = — , (5) 

we are left with the Hamiltonian 

H d /uj = tfb+^ + ^a z (tf + b)-^a x +5cr z + eV g , (6) 

where b and are the usual boson annihilation and cre- 
ation operators. This Hamiltonian can be solved nu- 
merically in a truncated boson Hilbert space, which is 
done below in order to determine the eigenenergies and 
the overlap factors between the empty and the occu- 
pied molecule. At a = 1 there is a cross-over to a 
regime where the fermion and boson degrees of freedom 
become strongly correlated due to the formation of po- 
laron states. When the levels are non-degenerate the po- 
laron leads to an increased localization of the electronic 
wavefunction. In both cases there is a strong reduc- 
tion of the effective tunneling amplitude between the two 
sites, which has consequences for conductance through 
the molecule. 



III. SEMI-CLASSICAL ANALYSIS 

In order to understand the polaron formation, we start 
by a semiclassical analysis of the spin-boson model. This 



is done in order to illustrate the physics and will not be 
used in the actual description described in the following 
sections. The semiclassical treatment is formally valid 
when uj <C t e ff ? but in the examples taken below this is 
not the case, and we instead resort to numerical diago- 
nalization of the Hamiltonian in Eq. (JJjJ). 

Treating the harmonic oscillator classically, the model 
can be solved in the electronic Hilbert-space. The elec- 
tronic Hamiltonian then becomes, 

fV^9X + 5 -2g/a \ 
H De =u>\ -2g/a -J2jjX-5\, (7) 

where X = x/i and I = yj\jmuj. For each of the two 
electronic eigenstates we thus have an effective poten- 
tial for the oscillator degree of freedom. These Born- 
Oppenheimer surfaces are 



Y± = l x 2 ± d2gX* + %+2y/TgX5 + P. (8) 
uj z V or 

In the electronic groundstate V_ there is only one mini- 
mum for a < 1 , whereas for a > 1 there are two minima. 
Thus a bifurcation occurs at a = 1. In the bifurcated 
domain, the electronic and bosonic degrees of freedom 
become highly correlated, because when the oscillator is 
localized in one of the two minima, the electron wavefunc- 
tion is changed accordingly. Thus the physics is similar 
to that of a small polarons (see Fig. EJ). In Fig. we 
show the bifurcated potential, the exact numerical eigen- 
values, and the effective tunneling coupling between the 
polaron states. From our exact diagonalization we thus 
find that the difference in energy between the ground- 
state and the first excited state becomes exponentially 
small in the polaron regime. 

In the semi-classical picture, the splitting of the po- 
laron states occurs because of tunneling. The tunneling 
amplitude is small, because in order to move, it must drag 
the oscillator displacement with it. Furthermore, a cou- 
pling to a dissipative environment will tend to localize the 
polaron even more, because it couples to displacement co- 
ordinate and thus tends to destroy the coherence between 
the two states, 19 resulting in a small effective tunneling 
coupling t e ff. When the effective coupling becomes very 
small, the two lowest eigenstates |^o) and form two 
polaron states which are approximate eigenstates: 

l^a P (M> = (IV>o) ± |V-i»/V2. (9) 

As long as the splitting in energy is larger than the tun- 
neling broadening of the levels, the polaron formation is 
not important for the transport other than it influences 
the Franck-Condon factors and the molecule can still be 
regarded as one quantum system (see Fig.|HJ). This situ- 
ation is analyzed in Section HVl However, if the effective 
tunneling rate t e g is smaller than the tunneling rates to 
and from the leads, a different physical picture emerges, 
because the dimer will behave as a "double dot" with 
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FIG. 2: Left panel: The semiclassical potential for g — 5 
and a = 2 and the exact eigenvalues (horizontal lines) found 
by numerical diagonalization of Eq. @. Note that the two 
lowest eigenvalues are almost degenerate, corresponding to 
two polaron states split by a small effective tunneling cou- 
pling t e ff. The effective coupling is shown in the right panel 
as a function of a together with the electron population of 
level a in the polaron state ^ ap , defined in Eq. @. The re- 
maining eigenenergies are not degenerate and therefore the 
corresponding eigenstates are delocalized on the molecule. 



weak interdot coupling but strong Coulomb interaction 
(see Fig. 0}. Thus a master equation treatment should 
take into account that the molecular system has more 
internal states, which is done in Section This regime 
is where the rectification occurs as discussed in detail 
below. 



IV. TRANSPORT PROPERTIES FOR T < t eff 

As argued, the transport properties of the dimer tran- 
sistor depend on the ratio between tunneling rates and 
internal time scales of the molecule. Let us first assume 
that all tunneling rates Y Ra ,T Rb ,T La and T Lb are much 
smaller than t e ff- Furthermore, throughout we consider 
the weak tunneling limit in the sense that the vibrational 
degrees of freedom are assumed to relax between tunnel- 
ing events, i.e. uo/Q ^> T, where Q is the quality factor 
of the vibrational mode. 11 If this limit is not satisfied the 
non-equilibrium vibron distribution should also be deter- 
mined, however, for the physics discussed here this is not 
important. Also, the broadening of the vibron sidebands 
is not included and assumed dominated by the thermal 
smearing such that k^T ^> uo/Q. To describe this situa- 
tion, we now use the master equation approach and only 
the occupations of the different charge states need to be 
determined by the kinetics. Allowing only two charge 
states, with occupations Po and Pi, the master equation 
is 



-Ar i - Wo = o, 

together with the normalization condition 
Po + Pi = 1. 



(10) 



(11) 



Here is the tunneling rate from state j to state i. The 
expression for the current then becomes 



J=(-e) 
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(r < t eff ). (12) 



In Eqs. (|l()|) . (|T2j) the tunneling rates are calculated us- 
ing Fermi's golden rule. When calculating the rates in 
Eq. Jl2j) we encounter the tunneling densities of states 
of the form: T 11 = 2tt \T kVij | 2 £(W, (v = L,R). 
The evaluation of the cross terms J2k T krj,a T kri,b^(^kr]) 
requires detailed knowledge of the relative phases of the 
tunneling amplitudes to the two parts of the dimer. For- 
tunately, a simplification is possible because Tj* a T k b ~ 
exp(z&<i), where d is the distance between the two parts 
of the dimer and k is the wave number of the lead elec- 
trons. Typically, one has kd ^> 1 and the cross-terms 
average to zero when summing over k. Thus, no inter- 
ference effects between the two tunneling paths onto the 
molecule is expected. With these observations, we obtain 



j=a,b i Q J 1 



2 e -PE i( 



and 



n v (E h -E iQ ), 
(13) 



,2 e 



a= E^'EK/oWi 

j=a,b i 1 J Q 

x [1 -n^- £/„)]■ 



(14) 



Here, |zq), (|/o)> denote initial (final) states of 
the empty (0) and filled (1) dimer, and 

n v (E) = l/(exp((£ - eV v )/k B T) + 1) (15) 

are the Fermi functions in the lead rj with applied voltages 
Vrj. Finally, E{ n and Ef n are the total dimer energies for 
the initial and final states with n electrons and Zq and 
Zi the corresponding partition functions. 

Because we are interested in the asymmetric situa- 
tion, we study the situation where one electrode cou- 
ples stronger than the other to the molecule, T R <C T L . 
Furthermore, we also allow for a skewed configuration 
as in Fig. ^ an d in this geometry, we have T Ra < T Rb 
and T La > T Lb . In Fig. we show examples of the 
differential conductance in the V — V g plane in the situa- 
tion where the intramolecular coupling between polaron 
states is strong, i.e. when T <C t e ff an d Eq. (|12|) applies. 
Furthermore, the bias voltage is applied symmetrically, 
so that V L -V R V/2. 



V. TRANSPORT PROPERTIES FOR T > t eS 

In the opposite case T ^> t e ff, the situation is very dif- 
ferent. In this case, an electron that tunnels onto level 
a does not resolve the tunneling splitting of the two po- 
laron states and the molecule no longer relaxes between 
tunneling and therefore the two polaron states can be 
considered as decoupled. In this case, we must treat the 
filled molecule as a quantum system having more possible 
states, since it can be occupied in one of the two polaron 
states, from now on denoted by a p and b p . Furthermore, 
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FIG. 3: (Color online) Left panel: Contour plot showing the 
differential conductance for the dimer molecule in the regime 
Eq. (|12f) where holds. The parameters are g — 5, a — 2, 



0, T* 



o, r Kb /r L 



0.04, r Lb /r L 



0.96, and 



T = O.lu. The bias voltage is applied symmetrically, so that 
Vl = —Vr — V/2. Notice the suppression of the differential 
conductance at low bias voltage which is due to a overlap 
between the filled and empty states. The right panel shows IV 
curves for the same parameters (but T — 0.05a;) for V g —V g o = 
(curve without gap at V = 0), —2oo (dashed) and 2ou. Here 
V g o defines the gate voltage at which the groundstates of the 
two charge states are degenerate. 



the molecule may be occupied in one of the delocalized 
states (see Fig. EJ. 

Moreover, we consider the following hierarchy of ener- 
gies ou > k^T > uj/Q > T > £ e ff- The last two inequali- 
ties imply that between tunnelings the molecule relaxes 
to a thermal distribution within each polaron subspace. 
Since the current drives the occupancies of the two po- 
laron states out of equilibrium, we need to determine 
the distribution functions of the two polaron states, P a 
and P 6p , respectively. If the molecule is occupied in a 
delocalized state it is in either polaron state with equal 
probability. The new set of master equations thus reads 
b P ) 



(with r = a p or r 



^r^oi ^0^10 — 0? 



— ^2 ^T^IO 



together with the condition 



P + £P r = l. 



(16a) 
(16b) 



(16c) 



The tunneling rates in Eq. ((16)) have contributions from 
both leads: 



E r "'' 

n=LM 



(16d) 



The tunneling rate for tunneling from lead r] into a p , is 
(again the oscillating cross-terms are ignored) 



i= a ' b i Q ,fi P 



2 e~^o 



— n v {E riP -E io ), 
(17) 



and similarly for tunneling into the polaron states b p . 
Here \fi P ) and \fi P ) mean the polaron states in Eq. ® 
if they are degenerate, i.e. (\fi~^) ± \f[ 2 ^)/y/2, where 
1/}^), i = 1,2 are the two degenerate states (meaning 
that eigenenergies differ by less than T). For the non- 
degenerate states, \fi P ) = \fi P ) = \fi) implying that 
the system ends up in either polaron state with equal 
probability. 

For the tunneling-out processes the rates are 



i= a > b h"J 

x [l-n v (E. ap -E fo )], 



(18) 



where the initial polaron states are defined in the same 
way as above. Similar expression holds for tunneling out 
of the polaron state b p . Thus, treating the dimer as sepa- 
rate polaron states, we find a different expression for the 
current: 



J=(-e) 
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(T > toff). 

(19) 



In Fig. 13 we show the results in the case where T ^> t e ^ 
and Eq. (|19|) holds. The IV characteristics are seen to 
be strongly asymmetric even though the bias is applies in 
a symmetric way, i.q.Vl — —Vr = V/2. The reason for 
the asymmetry is that the polaron state, which is weakly 
coupled to the right lead can block current through the 
other, more strongly coupled state. However, this block- 
ade only happens in one current direction, namely when 
the electrons run toward the weakly coupled electrode. 



VI. CROSS-OVER REGIME: GENERALIZED 
MASTER EQUATION 

In the regime where T and t e R are of the same order, 
we cannot use any of the two approaches above, because 
when two states in the molecular system are degenerate 
or closer than Y in energy, the current through the system 
has to be described as one coherent process. This is in 
general a difficult problem. However, in the limit of large 
bias voltage a generalized master equations2-i2ii22*2i has 
been developed. To study how the asymmetry of the 
large bias conductance, we here adopt this method. The 
large bias assumption implies that the voltage must be 
large compared to the temperature and the (effective) 
tunneling couplings as well molecular energies. The as- 
sumption of large bias leads to the a set of Markovian 
master equations for the density matrix Thus 
at large bias there is no correlation between individual 
tunneling events, but still the transfer of electrons is de- 
scribed coherently. The generalized master equation ap- 
proach can, in fact, describe the cross-over we discuss 
here, because the blocking effect found in Section V does 
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FIG. 4: (Color online) Left panel: Contour plots showing the 
differential conductance for the dimer molecule when Eq. (|19|) 
holds, for a = 2 (top) and a = 4 (bottom). Other parameters 
are g = 5, A = 0, F Ra = 0, T Rb /T La = 0.04, F Lb /F La = 0.96, 
and T = 0.1a;. The bias voltage is applied symmetrically, 
so that Vl — —Vr — V/2. Right panels show corresponding 
IV curves for different gate voltages as in Fig. |3| and with 
T — 0.05a;. Note the strong blocking of the current in one 
direction and the negative differential resistances. For a — 2 
the current in the blocked direction is larger than for a — 4, 
because the polaron is less localized on one electron level (see 

Fig. [g. 



not depend on time-correlations between tunnel events, 
but on how strongly a single electron is localized on the 
molecule during tunneling. 

Furthermore, in order to simplify the problem we con- 
sider only the two lowest states which correspond to the 
two polaron states located to the left or to the right with 
some effective coupling t e &. The generalized master equa- 
tions are then a modification of those developed in Ref.l2ll 
for a double dot system, because here we allow for in- and 
out tunnelings on both dots. For the case when the elec- 
tron current is from left to right (i.e. V < 0), we have 



dtpaa 


= T La P()Q _ T Ra Ppaa _ 


f- it e fi(Pba 


~ Pab), 


(20a) 
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(20d) 



The stationary solution for p is found by setting d t pij = 
and the current is then / = (—e)(p aa T Ra p + p bb T Rb p), 



which becomes for V < 

(-e)T L r R (x R T R + r^r^ /r R ) 



(r R + 2t l )x r t r + r^r^(i 



where 



^Ra v 



^Rb v 



^Lb v 



and where 



Xr 



4t e 



4A 2 + (r^) 2 



Y Rb v 

(21) 



(22) 



(23) 



contains all the t e ff and A dependencies. When Xr is 
large, the current reduces to the incoherent result of a 
dot with a doubly degenerate level: 



(-e)r L r^ 
+ 2r L 



(24) 



In the opposite limit of small Xr the time an electron 
spends on the dimer is too short for hopping between the 
sites to occur, and the current reduces to 



Iv 



<0 



(25) 



which is also the large V limit of Eq. ([TUj) or the sequen- 
tial tunneling limit of a double dot system in absence of 
interdot tunneling. 

The current for positive voltage is easily found from 
Eq. (|21[) by interchanging left and right, L <-> R. 

Let us now focus on the limit considered in the exam- 
ples in the previous two sections, namely (T R , t e ff) <^ T L 
and (r Ra ,t e ff) <C T^ 6 , and in these limits we find the 
following ratio between the currents for the two bias po- 
larizations 



Iv 



<0 



Iv 



>0 



Xr 



2 Xi 



(26) 



From this expression it is evident that degree of rectifica- 
tion is limited by the smallest tunneling-out rate through 
the right junction, T Rcip or the the inter-polaron tunnel- 
ing coupling, t e ff, through the parameter Xr m Eq. (|23|) . 
Eq. ((26)) also shows that when asymmetry is limited by 
Xr, i-e. when Y Rap /T R <C Xri increasing A in fact en- 
hances the rectification, because a finite A suppresses 
tunneling between the two sites. 



VII. DISCUSSION AND SUMMARY 

Lastly, we discuss the relevance for the experimental 
system in Ref. 0, where a single-electron transistor setup 
was made with single C140 molecules as the active ele- 
ment. Of course, our two level model cannot fully de- 
scribe the real molecule, nevertheless a comparison is 
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interesting. Signatures of a 11 meV intramolecular vi- 
brational mode were observed in Ref. |y and the coupling 
parameter g was estimated to be between 1 and 5jS4 If 
the polaron physics discussed here should have relevance 
for this or similar dimer systems, we must have a > 2, or 
correspondingly t ~ 50 meV, which is not unrealistic. In 
fact, rectification effects very similar to those predicted 
in the present paper have indeed been seen in some C140 
devices^ 

In conclusion, we have studied molecular system with 
an internal vibrational mode coupled to the charge dif- 
ference between two hybridized molecular levels giving 
rise to polaron formation. The polaron has distinct con- 
sequences for the transport properties, which we have 
calculated using different master equation approaches. 
First, we analyzed the case when the internal hybridiza- 
tion between polarons is larger than tunneling rates to 
the leads. This gives more or less symmetric IV char- 
acteristics, but with a number of vibron sidebands. Sec- 
ondly, we studied the case when the two polarons are 
weakly coupled and the tunneling rates to the leads dom- 
inate the kinetics. Is this case, another set of master 
equation was needed and highly asymmetric IV char- 



acteristics are predicted. Finally, we have analyzed the 
cross-over regime by a generalized master equation ca- 
pable of describing the coherent tunneling, but only for 
large bias voltages. 

The rectification mechanism suggested in this pa- 
per should be experimentally observable. Generally, 
the current- voltage characteristics of complex molecular 
transistors with strong coupling between charge and vi- 
brational degrees of freedom is a promising tool for study- 
ing the details of such devices. 
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